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Abstract 

The Higgs-boson production in association with a W-boson pair at e + e~ linear 
colliders is one of the important processes in probing the coupling between Higgs- 
boson and vector gauge bosons and discovering the signature of new physics. We de- 
scribe the impact of the complete electroweak(EW) radiative corrections of 0{a ew ) 
to this process in the standard model(SM) at the International Linear Collider (ILC), 
and investigate the dependence of the lowest-order(LO) and EW next-to-leading or- 
der(NLO) corrected cross sections on colliding energy yfs and Higg-boson mass. The 
LO and NLO EW corrected distributions of the invariant mass of W-boson pair and 
the transverse momenta of final W- and Higgs-boson are presented. Our numerical 
results show that the relative EW radiative correction^™) varies from —19.4% to 
0.2% when m H = 120 GeV and yfs goes up from 300 GeV to 1.2 TeV. 
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I. Introduction 



The Higgs-boson plays a very important role in the standard model(SM), it is responsible 
for the breaking of the electroweak(EW) symmetry and the generation of masses for 
the fundamental particles [H [2] . Unfortunately, it has not been directly detected yet in 
experiments. Searching for Higgs-boson within the SM and studying the phenomenology 
concerning its properties are the important tasks for the present and upcoming high 
energy colliders. Until now, we have obtained a lot of experimental data which agree 
with the predictions of the SM. These precise data constrain strongly the couplings of the 
gauge-boson to fermions and the couplings between vector gauge bosons, but present little 
information about the couplings between Higgs-boson and gauge bosons. Previous LEP 
experiments have provided the lower limit on the SM Higgs-boson mass as 114.4 GeV at 
the 95% confidence level, which is extracted from the results of searches for e + e~ — > Z°H° 
production^ H]. The present indirect evidences of the SM Higgs-boson mass through EW 
precision measurements indicate the 95% C.L. upper bound as m H < 182 GeV, when the 
lower limit on is used in determination of this upper limit [4J. 

Compared to hadron machine, e + e~ linear collider has the advantage of producing 
Higgs-boson in a particularly clean environment. That is essential for studying Higgs- 
boson properties since the cross section of Higgs-boson production is rather small due 
to the fact that the scalar Higgs-boson couples mainly to heavy particles. Furthermore, 
in the upper and lower limits for the SM Higgs-boson mass, H° decays mainly into the 
bottom-quark pair allowed by phase-space. Therefore, the properties of Higgs-boson with 
the mass in present limitation range, are very hard to be precisely probed at hadron 
colliders because of huge QCD backgrounds. 

The future International Linear Collider (ILC) is designed as a machine with the entire 
colliding energy in a range of 200 GeV < t/s < 500 GeV and an integrated luminosity 
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of around 500 in four years. The machine could be upgraded to y/s ~ 1 TeV 

with an integrated luminosity of 1 (a6) _1 in three years[5J. At this machine, Higgs- 
boson would be produced mainly via the Higgs-boson strahlung process, and the coupling 
of Higgs-boson to Z-bosons is probed best in the measurement of the cross sections of 
the Higgs-boson strahlung process e + e~ — > H°Z° and the WW/ZZ fusion processes 
e + e~ — > H°uu and e + e _ — > H°e + e~. The ILC would be an ideal machine to observe 
and study in detail an intermediate- mass or even a heavier Higgs-boson [SI El- It is sure 
that once the neutral Higgs-boson is discovered and its mass is determined, the HVV 
production through e + e~ — > H°VV (V = Z, W, 7) processes may provide the detail 
information of the coupling between Higgs-boson and vector gauge bosons, which directly 
reflects the role of the Higgs-boson in EW symmetry breaking. Moreover, a theoretical 
accurate estimation of these processes is essential, since e + e" — > H°VV (V = Z, W, 7) 
processes could be potential backgrounds for possible new physics. The investigation 
of the processes e + e _ — > H°W + W~ , e + e~ — > H°Z°Z° and e + e~ — > H°Z 0, y for testing 
the couplings between Higgs-boson and gauge bosons were studied at the tree-level in 
Ref . [9l [TO] . The coupling of guzz can be determined at a few percent level for a 120 GeV 
Higgs-boson with an integrated luminosity of 500 fb' 1 from the production cross section 
of the process e + e~ — > H°Z°Z° as shown in Ref. [8]. Recently, Y.-J. Zhou et al., studied 
the full 0(a ew ) EW corrections to the process e + e~ — > H°Z°Z°\X1\. They conclude that 
the 0(a ew ) EW radiative corrections to e + e~ H°Z°Z° process are in the range between 
— 15% and 1.0% when yfs varies from 400 GeV to 2 TeV, which should be taken into 
consideration in the future precise experiments. The process of Higgs-boson production 
in association with a pair of W-bosons at the ILC is another significant process. It is 
not only important in probing gnww coupling, but also possible to provide further tests 
for the quadrilinear couplings, such as HWWZ and HWW'j, which do not exist at tree- 
level in the SM, because these quadrilinear couplings would induce deviations from the 
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SM predicted observables. That means the precise theoretical prediction is necessary in 
experimental data analysis for e + e~ — ► H°W + W~ process. 

In this work we calculate the full one-loop EW corrections to the process e + e~ — > 
H°W + W~ in the SM. The paper is arranged as follows: In Section II we give the calcula- 
tion description of the Born cross section, and the calculation of full 0(a^ w ) EW radiative 
contribution to the e + e~ — > H°W + W~ process is provided in section III. In Section IV 
we present some numerical results and discussion, and finally a short summary is given. 

II. LO contributions of 0(a 3 ew ) to e+e~ -> H°W + W~ 
process 

In our LO and NLO calculations we apply FeynArts3.3 package [T2] to generate Feyn- 
man diagrams and convert them to corresponding amplitudes. The amplitude reductions 
are mainly implemented by applying FormCalc5.3 programs[13j. We use the 't Hooft- 
Feynman gauge in the leading-order calculation, if there is no other statement. The 
contribution to the process e + e~ — > H°W + W~ at the order of 0(a z ew ) in the SM is 
based on pure electroweak interactions. We ignore the contribution from the Feynman 
diagrams which involve electron-Higgs Yukawa coupling, because the concerning Yukawa 
coupling strength is proportional to the electron mass, and the electron mass is negligi- 
ble comparing with the colliding energy yfs at the ILC, Then there are eleven Feynman 
diagrams which we consider in the calculation for the process e + e~ — > H°W + W~ at the 
tree-level(shown in Fig{T]). The Feynman diagrams in Figfj] topologically can be divided 
into s-channel diagrams with intermediate 7- or Z°-boson, and t-channel diagrams with 
z/ e -exchanging. The notations for the process e + e~ — > H°W + W~ are defined as 

e + ( P i) + e-(p 2 ) - H°(p 3 ) + W + (p,) + W-(p 5 ), (2.1) 
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Figure 1: The full set of tree-level Feynman diagrams for the e + e — > H°W + W process 
neglecting the electron-Higgs coupling. 

where pi (i = 1 — 5) label the four-momenta of incoming positron, electron and outgoing 
particles, respectively. The differential cross section for the process e + e~ — > H°W + W~ 
at the tree-level is then obtained as 

dOtree = 7 X! \M tr ee\ 2 3 , (2.2) 
spin 

where Aitree is the amplitude of all the tree- level diagrams shown in Fig{TJ The factor | 
is due to taking average over the spins of the initial particles. d<&3 is the three-particle 
phase space element defined as 

(2,) 
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III. Virtual and real photon emission corrections to 

e + e~ — > H°W + W~ process 

The one-loop Feynman diagrams and relevant counterterm diagrams of process e + e~ — > 
H°W + W~ are also generated by using FeynArts3.3 package. At the EW NLO, there are 
as many as 2770 one-loop Feynman diagrams being taken into account in our calculation, 
and they can be classified into self-energy, triangle (3-point), box(4-point), pentagons- 
point) and counterterm diagrams. There are 82 pentagon diagrams. A representative set 
of pentagon diagrams can be found in FigfSJ The calculation of the one-loop diagrams has 
been performed in the conventional 't Hooft-Feynman gauge. The virtual contribution of 
C(«g W ) to e + e~ — > H°W + W~ process can be expressed as[ 



Acr v i rt ual = OtaeAirtual = , J, r- / ~d$ 3 ^ Re (Octree A4 virtual ) , (3.1) 

' ' spin 

where pi is the c.m.s. spatial momentum of the incoming positron. M. virtual represents the 
amplitude of 0(a^. w ) order Feynman diagrams including EW one-loop and counterterm 
diagrams. We use the dimensional regularization(DR) method to isolate the ultravio- 
let (UV) and infrared(IR) singularities. The numerical calculation of the one-, two-, three- 
and four-point integral functions are implemented by using the expressions presented in 
Refs.[15, 16J. The implementations of the scalar, vector and tensor five-point integrals are 
done exactly by using the approach presented in Refs. [TTf [TB] . The 0(a^ w ) calculation of 
the cross section for hard photon radiation process e + e~ — > H°W + W~ , y is accomplished 
by using Grace2.2.1 package [TU] . 

The total unrenormalized amplitude corresponding to all the one-loop Feynman di- 
agrams contains both ultraviolet (UV) and infrared (IR) divergences. The UV and IR 
divergences in loop integrals are isolated by adopting the dimensional regularization(DR) 
scheme. The relevant fields are renormalized by taking the on-mass-shell (OMS) scheme 
j. As we expect, the UV divergence contributed by virtual one-loop diagrams can 
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(1) (2) (3) (4) 




Figure 2: Some representative pentagon Feynman diagrams for the process e + e 
H°W + W- . 
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be cancelled by that contributed from the counterterms exactly both analytically and 
numerically 

In our one-loop calculation there exists soft IR divergence, but no collinear IR sin- 
gularity, because we keep nonzero electron(positron) mass. The soft IR divergence in 
virtual photonic corrections for the process e + e~ — > H°W + W~ can be exactly cancelled 
by adding the real photonic bremsstrahlung corrections to this process in the soft photon 
limit. In the real photon emission process, denoted as 

e+( Pl ) + e-(p 2 ) - H°(p 3 ) + W + (p 4 ) + W~(p 5 ) + 7 (p 6 ), (3.2) 

a real photon radiates from the electron/positron or W-bosons, and can be soft or hard. 
We adopt the general phase-space-slicing (PSS) method [2TJ to isolate the soft photon 
emission singularity part in the real photon emission process e + e~ — *> H°W + W~ , y . The 
bremsstrahlung phase space is divided into singular and non-singular regions, and the 
cross section of the e + e~ — > H°W + W~^ process is decomposed into soft and hard terms, 

A(T rea l = A(J so f t + Adhard = Vtreei^soft + Shard)- (3.3) 

We can consider that the radiated photon in the reaction e + e~ —>■ H°W + W~ , ~f can be ei- 
ther soft or hard, depending on the photon energy (Eg) in the center of mass system(c.m.s.) 
frame. The criterion is like that: if Eg < AE, the radiated photon is soft, otherwise it 
is hard, where AE is defined as AE = 5 s Et>, and Ef, is the electron beam energy in the 
c.m.s. frame and equals to y/s/2. Then theoretically both Aa so f t and Aahard should 
depend on the arbitrary soft cutoff 5 S , but the total EW one-loop correction(A<j to i) and 
Aa rea i should be cutoff 5 S independent. If the IR singularity in the soft photon emission 
process is cancelled exactly with that from the virtual photonic corrections, the indepen- 
dence of Aa to ti= Aa virtua i + Acr rea i) on the cutoff 5 S and fictitious photon mass m 7 , could 
be demonstrated in our calculculation. Since generally we take the soft cutoff 5 S to be a 
small value during our calculations, the terms of order AE / E^ can be neglected and the 
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soft contribution can be evaluated by using the soft photon approximation analytically 

[mug 



dAo- soft = -da tre 



«ew f d 3 p 6 



27T 2 J\p 6 \<AE 2p° 6 



Pi P2 P4_ P5 



(3.4) 



_Pi-Pe P2-P6 Pi-Pe P5-P6. 
where we define the four-momentum of radiated photon as = (p^po). As shown in 
Eg. (13.41) . the soft contribution has an IR singularity at AE = 0, which can be cancelled 
exactly with that from the virtual photonic corrections. The hard contribution is UV and 
IR finite, and can be computed directly by using the Monte Carlo technique. Finally, the 
corrected total cross section(<7 to< ) up to the order of 0(a^ w ) for the e + e~ — > H°W + W~ 
process is obtained by summing the 0{a^ w ) Born cross section(<Tt ree ), the 0{a^ w ) virtual 
cross section(A<7 v i rtua i), and the 0(af w ) cross section of the real photon emission process 
e+ e - ^ H°W + W~ 1 (Aa real ). 

Ctot = 0"trcc + A(Tt t = (Ttrcc + A(J v i rtua i + A(T rea \ = (Ttree (1 + Sew) j (^-5) 

where 5 ew = Virtual + ^soft + (Wd is the full 0(a ew ) electroweak relative correction. 



IV. Numerical results and discussion 

For the numerically verification of the gauge invariance in our LO calculation, we used 
Grace2.2.1 and FeynArts3.3/FormCalc5.3 packages separately, and adopted 't Hooft- 
Feynman gauge and unitary gauge respectively in calculating the Born cross section of the 
process e + e~ — > H°W + W~ , and got the coincident numerical results. We use also both 
Grace2.2.1 and FeynArts3.3/FormCalc5.3 packages to calculate the tree-level cross section 
of process e + e~ — > H°W + W~ in 't Hooft-Feynman gauge with the input parameters used 
in Ref.JTD], and compare our numerical results with those in Ref.JTD] in order to verify 
again the correctness of our tree-level calculation. The input parameters used in Ref. [TU] 
are taken as: m z = 91.18 GeV, mw = 80.1 GeV, a ew {m z ) = 1/128, y/s = 500 GeV and 
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m H {GeV) 


(T tree (/6)(Ref.[TDl) 


cr tree (/6) (Grace) 


o~tree(fb) (FeynArts) 


120 


5.63 


5.6378(6) 


5.6374(6) 


150 


3.51 


3.5123(3) 


3.5126(3) 


200 


1.47 


1.4714(1) 


1.4715(1) 



Table 1: The comparison of the numerical results of the tree-level cross sections 
for the process e + e~ — > H°W + W~ with the corresponding ones in Ref.[10j, by 
taking the same input parameters as in Ref.fTU] and using Grace2.2.1 and Fey- 
nArts3.3/FormCalc5.3 packages separately. 

an effective sm 2 9w = 0.232. The numerical results of the tree-level cross section for the 
process e + e~ — > H°W + W~ are listed in Table [TJ We can see that Table [T] yields results 
in mutual agreement between ours and those presented in Ref.|10j. 

The following numerical computation is performed in the a ew -scheme and the initial 
EW input physical parameters are taken as [14] : 

m w = 80.403 GeV, m z = 91.1876 GeV, a ew {0) = 1/137.03599911. (4.1) 

The charged leptons have the mass values as 

m e = 0.51099892 MeV, = 105.658369 MeV, m T = 1776.99 MeV. (4.2) 

For the quark masses, beside the top-quark mass m t = 172.5 GeV, we take the set 
m u = md = 66 MeV, m s = 150 MeV, m c = 1.5 GeV, and = 4.7 GeV. Since 
quark-mixing effects are suppressed, we set the CKM matrix to the unit matrix. 

Except the input parameters shown above, we have to fix the values of the IR regulator 
m 7 , the fictitious photon mass, and soft cutoff 5 S = AE/E^ during our one-loop numerical 
calculation. In fact, if the one- loop calculation is correct and the IR divergency is really 
cancelled, the total cross section should be independent of these two parameters. Our 
numerical results show that the cross section contribution at 0(aj w ) order, Aa tot = 
Acr re ai + Aa v i rt uai, * s invariant within the calculation errors when the fictitious photon 
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Figure 3: (a) The 0(a^ w ) contribution parts of cross section for e + e~ — > H°W + W~ 
process as the functions of the soft cutoff 5 S = AE/Ej, in conditions of m 7 = 10~ 10 GeV, 
ran = 120 GeV and y/s = 500 GeV. (b) The amplified curve for Aa tot of Fig.3(a) versus 
S s including calculation errors. 

mass m 7 varies from 10~ 15 GeV to 10 _1 GeV in conditions of 5 S = 10 -3 , vtlh = 120 GeV 
and y/s = 500 GeV. 

As a check of the correctness in our calculation of radiative correction, we present 
Figs(3](a-b) to show the independence of the total 0(a^ w ) electroweak contribution to 
e + e~ — > H°W + W~ process on soft cutoff 6 S , assuming m 7 = 10 -5 GeV, mu = 120 GeV 
and y/s = 500 GeV. As shown in Figj3]^a), both Acr v i rtua i + Acr so f t and AcThard strongly 
depend on soft cutoff S s , but the total 0(a^ w ) EW relative radiative contribution A<r tot = 
A<T V i r tuai+A<7 re ai is cutoff 5 S independent within the range of calculation errors as expected. 
In order to show the curve of Aatot more clearly, we present the amplified curve for Acr toi 
including calculation errors in Figl^b). In further calculations, we fix m 7 = 10~ 5 GeV 
and 5 S = 10~ 3 . 

In FigJU^a) we present the curves of the LO and one-loop EW corrected cross sections as 
the functions of colliding energy y/s, with uih = 120 GeV, 150 GeV, 200 GeV, separately. 
FigJU^b) shows their corresponding relative EW radiative corrections(5 e ^ = ^°" iot ) for 
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Figure 4: (a) The LO and EW one-loop corrected cross sections for the process e + e~ — > 
H°W + W~ as the functions of colliding energy y/s with m# = 120 GeV^, 150 GeV, 
180 GeK respectively, (b) The corresponding relative EW radiative corrections (S ew ) 
versus y/s. 

the data presented in FigJH(a). We find from FigsJUa-b) that the LO cross sections 
are sensitive to the Higgs-boson mass and generally suppressed by the one-loop EW 
corrections except for ran = 120 GeV and the colliding energy being in the range of y/s > 
630 GeV. FigHJb) shows that the relative radiative correction due to full EW one-loop 
contributions in the y/s vicinity close to the threshold of H°W + W~ production, becomes 
rather large. That is because of the Coulomb singularity effect from the diagrams involving 
the instantaneous virtual photon exchange in loop which has a small spatial momentum. 
To show the numerical results as presented in Figs|H(a-b) more precisely, we list some 
typical numerical results of the tree-level, one-loop EW corrected cross sections and the 
relative EW radiative correction(5 eu , = Aa to t I ®tree) for the process e + e~ — > H°W + W~ in 
Table [2J There they are in conditions of <Js = 400 GeV, 600 GeV, 800 GeV, 1000 GeV, 
1200 GeV and m H = 120 GeV, 150 GeV, 180 GeV separately. We can read out also 
from Fig 11(b) that when m# = 120 GeV and y/s goes up from 300 GeV to 1.2 TeV, the 
relative EW radiative correction, S ew , increases from —19.4% to 0.2%. 

In Figj5|a) we present the plot of the LO and EW one-loop corrected cross sections 
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(GeV) 


m H (GeV) 


CTtree(fb) 


vtot(fb) 


M%) 




120 


4.1676(4) 


3.953(2) 


-5.15(5) 


400 


150 


1.9144(2) 


1.719(1) 


-10.21(5) 




180 


0.70438(7) 


0.5823(5) 


-17.33(7) 




120 


5.0945(5) 


5.085(3) 


-0.19(6) 


600 


150 


3.5357(3) 


3.392(2) 


-4.06(6) 




180 


2.4562(3) 


2.254(2) 


-8.23(8) 




120 


4.3205(6) 


4.354(3) 


0.78(7) 


800 


150 


3.3009(4) 


3.202(2) 


-3.00(6) 




180 


2.5672(3) 


2.393(2) 


-6.79(8) 




120 


3.5817(6) 


3.607(3) 


0.71(8) 


1000 


150 


2.8635(4) 


2.776(2) 


-3.06(7) 




180 


2.3418(3) 


2.182(2) 


-6.82(8) 




120 


3.0044(5) 


3.010(3) 


0.19(8) 


1200 


150 


2.4678(4) 


2.380(2) 


-3.56(8) 




180 


2.0770(3) 


1.926(2) 


-7.27(9) 



Table 2: The numerical results of the tree-level, one-loop EW corrected cross sections 
and the relative EW radiative correction(<5 ew = Aa tot / 'cr tree ) for the process e + e" — > 
H°W + W~ , by taking = 400 GeV, 600 GeV, 800 GeV, 1000 GeV, 1200 GeV 
and m H = 120 GeV, 150 GeV, 180 GeV separately. 
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m (GeV) 



Figure 5: (a) The LO and EW one-loop corrected cross sections for the process e + e~ — > 
H°W + W~ as the functions of the Higgs-boson mass m# with \/s = 500 GeV, 1000 GeV. 
(b) The corresponding relative EW radiative corrections versus Higgs-boson mass. 

for the process e + e~ — > H°W + W~ as the functions of Higgs-boson mass uih by taking 
\fs = 500 GeV and 1000 GeV separately, and their corresponding relative EW radiative 
corrections versus Higgs-boson mass are depicted in FigJI^b). From these two figures 
we can see again the cross sections at LO and up to one-loop order for the process 
e + e~ — > H°W + W~ are all sensitive to the Higgs-boson mass as already shown in Figs|H(a- 
b). And both curves for y/s = 500 GeV and 1000 GeV in FigJS^b) demonstrate that the 
absolute relative EW radiative corrections for yfs = 500 GeV and 1000 GeV go up quickly, 
when Higgs-boson mass goes up from from 100 GeV to 180 GeV, but become relative 
stable when m# runs from 180 GeV to 200 GeV. 

Since we consider the CP-conserving SM, the distributions of transverse momenta of 
VF~-boson should be the same as that of . Therefore, we shall not provide the distri- 
bution of but only the . We depict the differential cross sections of transverse 
momentum of if°-boson at the LO and up to NLO (daio,NLo/dPT ) m FigJ51( a )! an d the 
distributions of daLo/dp^ + and da^Lo/ 'dp^ + in FigJ^b) separately, in the conditions of 
ttlh = 120 GeV and y/s = 500 GeV. We can see from Figsj5]Ja-b) that when the trans- 
verse momentum p^°{p^ + ) is in the range smaller than 90 GeV, the LO differential cross 
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50 100 150 200 50 100 150 200 



p"(GeV) p"(GeV) 

Figure 6: The distributions of the transverse momenta of H°- and VT + -boson for the 
e + e~ — > H°W + W~ process at the LO and up to EW one-loop level with yfs = 500 GeV 
and run = 120 GeV. (a) for i7°-boson, (b) for iy + -boson. 

section of daLo/dPr {d (T Lo/dp^ + ) is slightly enhanced by the EW one-loop corrections, 
while it is suppressed when Pt°(Pt + ) is larger than 100 GeV. 

We plot the distributions of W-pair invariant mass, denoted as Mww, & t the LO and 
up to EW one-loop level in FigJ7|by taking tuh = 120 GeV and \/s = 500 GeV. We can 
see from the figure that when My/w < 270 GeV the one-loop EW corrections enhance 
the LO differential cross section daio/dMww slightly, while in the Mww region being 
larger than 280 GeV the daio/dMww is obviously suppressed by the EW corrections. 

V. Summary 

The future e + e~ linear collider, International Linear Collider (ILC), will probably offer 
the cleanest environment to probe the SM more precisely and discover the signature of 
new physics. In this paper we have shown that the phenomenological effects due to the 
contribution from complete one-loop EW terms of 0(a^ w ), can be demonstrated in the 
study of the electron-positron into Higgs-boson in association with a W-boson pair in 
the standard model(SM) for all energies ranging from 300 GeV to 1.2 TeV at ILC. We 
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Figure 7: The distributions of the invariant mass of W-pair at the LO and up to EW 
one- loop level when m# = 120 GeV and y/s = 500 GeV. 

discuss the dependence of the effects coming from the EW NLO contribution to the cross 
section of process e + e~ — > H°W + W~ , on colliding energy y/s and Higgs-boson mass. 
We present the LO and EW NLO corrected distributions of the transverse momenta of 
final particles and the differential cross section of H^-pair invariant mass. We find that 
the LO cross section of process e + e _ — > H°W + W~ is sensitive to Higgs-boson mass and 
generally suppressed by the one-loop EW corrections except when m H = 120 GeV and 
y/s > 630 GeV. Our numerical results show that when m H = 120 GeV and y/s goes up 
from 300 GeV to 1.2 TeV the relative EW radiative correction varies from —19.4% to 
0.2%. It shows that the EW one-loop radiative corrections to e + e" — > H°W + W~ process 
can be significant and should be included in any reliable analysis. 
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